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DEVELOPMENT  OF  A  TARGET  EVASION  REGION 
HAVING  A  CONDITIONAL  GAUSSIAN  DISTRIBUTION 


I.  INTRODUCTION 


In  an  earlier  report  (reference  1),  a  generalized  probability  density 
function  was  developed  for  modeling  maneuvering  target  position  inside  a 
bounded,  time-varying,  evasion  region  (reference  2).  The  approach  taken  was 
to  assume  beta  density  functions  for  changes  in  evasion  speed  and  course,  and 
then  formulate  the  resultant  positional  density  function.  Beta  density 
functions  were  employed  because  of  their  flexibility;  i.e.,  uniform,  ramp,  or 
other  type  density  functions  can  be  easily  developed  via  proper  selection  of 
the  shaping  parameters. 

One  reason  for  developing  these  probabilistic  evasion  regions  was  to 
formulate  alternate  firing  strategies  for  advanced  weapons.  The  mean  or  mode 
of  the  density  function  can  serve  as  an  aimpoint  based  on  target  evasion. 
Another  use  for  these  density  functions  is  in  the  area  of  online  measure-of- 
ef f ectiveness  (MOE)  algorithms.  MOE  algorithms  allow  real-time  prediction  of 
weapon  performance  in  a  tactical  encounter.  To  be  of  value,  these  algorithms 
must  take  into  account  the  initial  target  motion  analysis  (TMA)  uncertainty 
region,  the  target  evasion  region,  weapon  delivery  errors,  and  weapon  destruct 
characteristics . 

The  work  presented  in  reference  1  provides  an  analytical  formulation  of 
the  evasion  statistics;  but  these  statistics  must  be  combined  with  the  TMA 
uncertainty  region,  which  has  Gaussian  characteristics.  This  results  in  a 
total  target  uncertainty/evasion  region  density  function.  The  combined  region 
is  very  difficult  to  determine  analytically  because  of  the  complexity  of  the 
evasion  region  density  functions.  However,  if  it  is  assumed  that  both  the  TMA 
uncertainty  region  and  the  evasion  region  are  Gaussian  distributed,  then  the 
two  densities  could  very  easily  be  combined  to  form  another  Gaussian  density 
function  (reference  3).  Moreover,  this  combined  density  function  would 
contain  all  of  the  pertinent  statistical  information  inherent  in  the  TMA 
uncertainty  and  the  evasion  region  densities. 

This  report  documents  development  of  a  Gaussian-distributed  target 
evasion  region.  To  obtain  such  a  region,  target  evasion  speed  is  modeled  by  a 
Rayleigh  distribution,  while  target  evasion  course  uses  a  uniform  density 
function.  Although  the  uniform  course/Rayleigh  speed  model  yields  a 
positional  density  function  that  is  Gaussian  distributed,  there  is  a 
disadvantage  in  using  the  Rayleigh  density  function  because  of  its  infinite 
tail.  Since  maximum  evasion  speed  (Smax)  is  finite,  the  modeling  parameter 
in  the  Rayleigh  density  that  determines  the  most  likely  evasion  speed  must  be 
less  than  Smax.  In  fact,  it  has  been  shown  (reference  4)  that  this  modeling 
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parameter  should  be  less  than  0.41  Smax  to  ensure  that  there  is  less  than  a 
5  percent  chance  of  obtaining  evasion  speeds  greater  than  To  overcome 

this  constraint,  a  conditional  Rayleigh  distribution  must  be  used.  The 
resulting  positional  uncertainty  region  then  becomes  a  conditional  Gaussian 
density  function.  Finally,  it  is  shown  that  a  beta  density  can  be  employed  to 
approximate  a  conditional  Rayleigh  distribution  for  large  values  of  the 
Rayleigh  modeling  parameter. 
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II.  EVADING  TARGET  POSITIONAL  DENSITY  FUNCTIONS 


Two  analytical  density  functions  that  model  the  evasion  region  of  a 
target  were  determined,  along  with  their  respective  means  and  variances.  The 
first  positional  density  was  formulated  using  uniform  course/Rayleigh  speed 
density  functions  as  inputs.  The  second  was  developed  for  a  uniform 
course/truncated  Rayleigh  speed  combination.  In  this  section,  the  analytical 
formulations  of  the  models  are  presented  and  evaluated. 


UNIFORM  COURSE/RAYLEIGH  SPEED  COMBINATION 

The  form  of  a  Rayleigh  density  function  to  model  evasion  speed  is  given 
by  (reference  5) 


fs(s)  -  5  expt-^m,2))  _  „  ,  s  < 

a 

where  S  is  the  target  speed.  The  modeling  parameter  a,  which  determines  the 
mode  of  the  density,  would  represent  the  most  likely  evasion  speed.  In 
addition,  this  parameter  appears  in  the  mean  value  of  the  density,  which  is 
written  as  (reference  5) 


a/rf 

mS  ~  S2  * 


(2) 


Figure  1  shows  a  Rayleigh  density  function  for  modeling  target  evasion 
speed.  The  mode  of  the  density  occurs  at  a,  which  was  set  to  a  value  of 
10.  Included  in  this  illustration  are  the  probabilities  of  obtaining  evasion 
speeds  greater  than  Smax  (maximum  target  speed)  for  Smax  values  of  25,  30, 
and  40  knots.  As  can  be  seen,  even  for  a  maximum  speed  of  25  knots,  the 
probability  of  exceeding  this  value  is  low.  But  a  represents  the  most 
likely  evasion  speed,  which  is  10  knots  in  this  example.  For  large  values  of 
a,  the  probability  of  obtaining  evasion  speeds  greater  than  Sma:iC 
increases.  This  probability  is  determined  from 


P(S  _>  Smax)  -  exp{-S^ax/ (2a^)} . 


(3) 


Since  the  target  evasion  regions  are  given  as  a  function  of  time,  changes  in 
target  speed  can  be  scaled  into  changes  in  distances  by 

r  =  st,  (4) 

where  t  is  a  deterministic  quantity  representing  evasion  time  (growth  of  the 
evasion  region  (reference  2)).  Now,  the  density  function  for  speed  changes 
can  be  transformed  into  an  evasion  range  density  function  of  the  form 
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Figure  1.  Rayleigh  Density  Function  (a  =  10) 


£  (r)  >  ,  o  <  r  < 

^  /si  ^  0 

0C  t 


and 


(5a) 


at/ir 

mr  =  ST  • 

For  the  evasion  course  9,  a  uniform  density  of  the  form 
f@(9)  =  l/2ir,  — ir  <  9  <  ir, 

is  used.  Since  fr(r)  and  f@(9)  are  independent,  the  joint  density  frQ(r,9) 
becomes 


(5b) 


(6) 


f  fr  _  r  exp{-r2/(2a2t2)} 

frQ(r,9)  -  22 

2ira  t 


0  <  r  < 


-ir  <  9  <  ir. 


(7) 
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To  transform  equation  (7)  from  polar  to  rectangular  coordinates,  frQ(r,0) 
is  multiplied  by  the  Jacobian  (reference  5),  which  is  given  by  (reference  1) 


J  =  1 +  y2  .  (8) 

Transforming  the  joint  density  function  into  rectangular . coordinates  yields 


(  )  _  exp{-(x2  +  y2)/(2a2t2)} 

VX’y)  -  2ffa2t2 

which  is  a  joint  Gaussian  probability  density  function  with  variance 


(9a) 


ax  =  °y  = 


(9b) 


Equation  (9a)  is  an  analytical  description  of  the  statistics  of  the  target 
position  based  on  the  evasion  time  t  and  modeling  parameter  a.  This  type  of 
evasion  strategy  can  be  easily  combined  with  the  initial  Gaussian  uncertainty 
region  obtained  from  TMA  techniques.  An  analytical  formulation  for  this 
combined  uncertainty  region  is  a  critical  element  in  any  MOE  algorithm  for 
realistic  quantification  of  weapon  performance. 

From  equation  (9b),  it  is  seen  that  the  variance  for  the  positional 
density  function  is  determined  from  the  parameter  at,  which  also  determines 
the  mode  and  mean  of  the  Rayleigh  density  given  by  equation  (5a).  As  was 
shown  in  figure  1,  a  must  be  small  to  minimize  the  probability  of  obtaining 
evasion  speeds  greater  than  Smax.  This  requirement  results  in  a  narrow 
positional  density  function  to  describe  the  evasion  region  (see  the 
three-dimensional  plot  in  figure  2).  As  a  increases,  the  variance  given  by 
equation  (9b)  increases,  causing  the  positional  density  function  to  spread 
out.  Figures  3,  4,  and  5  are  plots  of  the  positional  density  functions  for 
a  =  15,  20,  and  30,  respectively.  Also  contained  in  these  figures  is  the 
probability  of  obtaining  speeds  greater  than  SmaxN  As  can  be  seen,  even  for 
a  =  15,  there  is  a  high  probability  of  obtaining  evasion  speeds  greater  than 

^max* 


If  an  instantaneous  motion  model  (course  and  speed  changes  are  step 
functions)  is  used  to  update  evading  target  position,  then  the  evasion  region 
would  be  a  circle  centered  at  the  initial  alerted  target  position  with  a 
radius  given  by 


rmax 


^max^ • 


(10) 


For  large  values  of  a,  the  resultant  positional  density  function  for  this 
evasion  region  assigns  a  non-zero  probability  to  points  outside  the  circle  of 
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*(x,y) 


Figure  2.  Three-Dimensional  Gaussian  Positional  Density  Function  (a  =  10) 


radius  rmax.  This  becomes  apparent  by  examining  the  variance  of  the 
positional  density  function  given  by  equation  (9b).  The  three-sigma  region  is 
3at  and,  if  a  >>  Smax/3,  the  three-sigma  region  would  be  larger  than 
rmax.  Accordingly,  a  positional  probability  region  that  is  Gaussian 
distributed  can  only  model  targets  evading  at  slow  evasion  speeds.  Otherwise, 
such  a  model  does  not  accurately  depict  the  real  evasion  region. 
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Figure 


Figure 


4.  Three-Dimensional  Gaussian  Positional  Density  Function  (a  =  20) 


3.  Three-Dimensional  Gaussian  Positional  Density  Function  (a  =  15) 
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f(x,y) 
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Figure  5.  Three-Dimensional  Gaussian  Positional  Density  Function  (a  =  30) 
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CONDITIONAL  RAYLEIGH  DENSITY  FUNCTION 


Because  of  the  restriction  placed  on  the  evasion  speed  density  function 
in  the  preceding  section,  the  resultant  Gaussian  uncertainty  region  can  be 
used  only  for  targets  evading  at  slow  speeds.  To  overcome  this  limitation  and 
still  obtain  a  Gaussian  evasion  region,  a  truncated  or  conditional 
distribution  must  be  employed.  This  density  is  found  from  (reference  5) 


f o (s | s  <  s  )  =  — 

q  1  -  max  S 


fs(S) 


j^max 


fs(S)dS 


where  Smax  is  the  maximum  evasion  speed.  Employing  equation  (11),  the 
resultant  conditional  Rayleigh  density  function  becomes 


(ID 


Us  s  <  s  )  = 

S  1  -  max 


 S  exp{-S2/(2q2)} 


a2F(S  ) 
max 


(12) 


where 


F(smax)  =  1  ~  exp{-s2ax/(2a2)}. 


(13) 


Once  again,  using  equation  (4),  a  conditional  density  function  for 
evasion  range  can  be  obtained: 


f  (r | r  <  r  )  = 
r  1  -  max 


r  exp{-r2/(2q2t2)> 

a2t2F(S  ) 
max 


(14) 


Plots  of  the  truncated  Rayleigh  density  function  given  by  equation  (12)  are 
shown  in  figures  6a  through  6d  for  a  =  10,  15,  20,  and  30,  respectively.  In 
all  of  these  plots,  the  Rayleigh  density  was  truncated  at  Smax  =  35  knots. 

As  can  be  seen  from  the  plots,  for  small  values  of  a,  the  truncation  has 
negligible  effect  but,  as  a  increases,  more  of  the  tail  is  dropped,  causing 
the  height  near  the  tail  to  increase  to  keep  the  total  area  equal  to  one. 

Thus,  the  conditional  speed  density  function  places  more  weight  on  the  evasion 
speeds  between  zero  and  Smax. 
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CONDITIONAL  GAUSSIAN  DENSITY  FUNCTION 


Using  the  uniform  density  function  given  by  equation  (6)  to  model  the 
evasion  course,  the  joint  Gaussian  probability  density  function  for  the 
evading  target  position  becomes 


f 

xy 


(x.yjx" 


2 

+  y 


< 


2 

r 

max 


)  = 


exp{-(x2  +  y2)/(2a2t2)} 

2ira2t2  F(S  ) 
max 


(15) 


Thus,  the  positional  density  function  for  the  target  evasion  region,  equation 
(15),  is  a  conditional  Gaussian  probability  function,  where  the  normalization 
factor  is  F(Smax).  The  error  function  that  would  normally  be  associated 
with  a  conditional  Gaussian  density  function  does  not  appear  because  the 
bounding  region  is  a  circle  of  radius  rmax.  Therefore, 


F(r  )  = 
max 


If. 


2  2.. 

liraz  2  2  2 

x  +y  <r 

-  max 


exp{-(x2  +  y2)/ (2a2t2)}dxdy , 


(16) 


which  reduces  to 


F(r  )  »  1  -  exp{-r2  /(2a2t2)}  =  F(S  ), 
max  r  max  max 


(17) 


since  rm^y  —  Smaxt. 


The  mean  value  of  equation  (15)  is  zero,  and  the  variance  is  given  by 


2 

a  = 
x 


f f  x2fxy(x’ylx2  +  y2  1  rmax)dxdy ’ 


xy 


2  2  2 
x  +y  <_r 


max 


(18) 


2 

a  = 

y 


ff  y2fxv(x’yl*2  +  y2 <  rLx)dxdy- 


xy 


2  2  2 
x  +y  <r 

-  max 


To  determine  equation  (15)  is  substituted  into  equation  (18)  and  a 
change  of  variables  to  polar  coordinates  is  made,  yielding 


(19) 


r  2tt 
max 


2 

ax  -  2  2 

2irct  t  F(S  )  n 
max  0 


- - f  f  r^sin^©  exp{-r2/(2a2t2)}d0dr . 


(20) 
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Integrating  equation  (20)  with  respect  to  9  results  in 


max 


2 

a  = 


X  2a2t2F(S  ) 

max  0 


- -  j  r^  exp{-r^/(2a^t^)}dr . 


(21) 


By  making  the  substitution  z  =  r^,  equation  (21)  becomes 

2 


max 


<7  = 


X  4a2t2F(S  )  n 
max  0 


- -  f  z  exp{“Z^/(2a^t2)}dz . 

Wc  J 


(22) 


As  shown  in  reference  5,  equation  (22)  reduces  to 


2  2  2  2 
o  0  0  r  exp{-r  /(2a  t  )} 
2  2  2  max  max  J 

%  *  a  11 - 2HS - ) - 

max 


(23) 


Similarly,  cr^  can  be  evaluated,  resulting  in 


9  „  .  r2  exp{-r2  /(2a2t2)} 

2  2^  2  max  r  max  J 

°y  -  *  c  -  - ms — > - 

max 


(24) 


Comparing  equation  (9b)  with  (23)  and  (24),  it  is  seen  that  the  difference 
between  the  variance  of  the  two  positional  density  functions  —  equations  (9a) 
and  (15)  —  is  in  the  second  term.  Recall  that  the  conditional  probability 
density  function  was  developed  to  keep  the  maximum  evasion  speed  finite 
(Smax  <  °°)*  So,  ^  smax  *  equations  (23)  and  (24)  reduce  to 
equation  (9b);  i.e.,  the  second  term  goes  to  zero.  To  show  this,  let 


K  =  lim 
S 

max 


rr2  exp{-r2  /(2a2t2)} 
max  r  v  max _ _ 

2F(S  ) 
max 


which  can  be  written  as 


K  =  lim 

S  *°° 
max 


s2  t2 

max 


2F(S  '  )  exp{S2  /(2a2)} 
max  r  w  max  J 


(25) 


(26) 


Substituting  equation  (13)  into  equation  (26)  and  employing  L'Hospital's  rule 
results  in 


13 


=  0. 


(27 


K  =  lim 

S  -*30 
max 


2  2 
2a  t 


exp{Smax/(2a2)>. 


Figures  7a  through  7d  show  plots  of  the  conditional  positional  density 
function  for  the  values  of  a  given  in  figures  6a  through  6d.  As  can  be  seen 
from  these  plots,  for  large  values  of  a  the  density  functions  flatten  out 
since  a  is  directly  proportional  to  the  variance  (see  equations  (23)  and 
(24)).  Also,  the  positional  density  function  goes  to  zero  for  r  _>  rmax. 

Thus,  the  goal  of  obtaining  a  Gaussian  density  function  to  model  the 
evasion  region  and  still  not  have  restrictions  on  the  evasion  speed  density 
function  has  been  achieved. 
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Figure  7a.  Conditional  Gaussian  Density  Function  (a  =  10) 
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APPROXIMATION  TO  A  RAMP  DENSITY  FUNCTION 


Figure  6d  shows  that  as  a  increases,  the  truncated  Rayleigh  density 
function  begins  to  resemble  the  ramp  density  function  (reference  1).  In  fact, 
it  can  be  shown  that  as  a  +  00 ,  the  truncated  Rayleigh  density  function 
becomes  a  ramp  density  function  and,  from  reference  1,  a  ramp  speed  density 
function  combined  with  a  uniform  course  density  function  yields  a  uniform 
positional  density  function.  To  show  that  the  truncated  Rayleigh  function 
approximates  a  ramp  function  as  a  -»  equation  (12)  can  be  written  as 


fS<sls 


<  S  )  = 
-  max 


2 

a 


S  exp{-S2/(2a2)} 

[1  -  exp{-S2  /(2a2)}] 


Expanding  the  exponentials  into  a  series  yields 


(28) 


f„(s|s  <  S  )  =  - 

5  1  -  max  *" 


1  - 


2a 


4a 


max 


s4  s6 

max  max 
— —  +  — ; — 


(29) 


-  2  4a  6a 
Now  take  the  limit  as  a  *>  °°,  and  equation  (29)  reduces  to 


fc(S  S  <  S  ) 
S  -  max 


2S 

S 

max 


(30) 


Figures  8a  through  8d  are  plots  of  equation  (12)  for  a  =  50,  100,  200, 
and  300,  respectively;  figure  9  is  a  plot  of  a  ramp  density  function.  A 
comparison  of  figures  8  and  9  shows  that,  even  for  a  =  100,  the  truncated 
Rayleigh  function  is  a  good  approximation  of  the  ramp  density  function. 


Figures  10a  through  lOd  are  plots  of  the  positional  density  functions  for 
the  values  of  a  in  figures  8a  through  8d.  These  plots  show  that  the  center 
flattens  out  as  a  increases  and  finally  approximates  a  uniform  positional 
density  function  at  a  =  300. 
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Figure  8a.  Truncated  Rayleigh  Density  Function  (a  =  50) 


18 


19 


20 
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Figure  10a.  Truncated  Gaussian  Density  Function  (a  =  50) 
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Figure  10b.  Truncated  Gaussian  Density  Function  (a  =  100) 
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Figure  10c.  Truncated  Gaussian  Density  Function  (a  =  200) 

4  f(x,y) 


Figure  lOd.  Truncated  Gaussian  Density  Function  (a  =  300) 
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III.  CONCLUSIONS 


This  report  has  examined  various  probability  density  functions  for 
modeling  target  evasion  regions.  Previous  work  (reference  1)  employed  beta 
density  functions  exclusively  to  achieve  this  goal.  The  work  presented  here 
uses  a  Rayleigh  density  function  to  model  speed  changes  and  a  uniform  density 
function  to  model  evasion  course.  Unlike  the  beta  density  function  with  its 
finite  limits,  the  Rayleigh  density  function  has  an  infinite  tail.  Therefore, 
to  obtain  evasion  speeds  that  are  finite,  the  Rayleigh  density  function  was 
truncated,  resulting  in  a  conditional  probability  density  function. 

The  attractiveness  of  the  Rayleigh  speed/uniform  course  combination  is 
that  the  concomitant  positional  density  function  for  the  target  evasion  region 
is  Gaussian.  The  Rayleigh  distribution  has  a  parameter  (a)  that  is 
analogous  to  the  beta  density  function's  shaping  parameters;  a  determines 
the  mode  of  the  density,  as  well  as  its  mean  and  standard  deviation.  This 
parameter  also  appears  in  the  variance  of  the  resulting  Gaussian  positional 
distribution.  Consequently,  the  spread  of  the  Gaussian  distribution  is 
directly  related  to  a.  The  results  obtained  in  this  study  also  show  that  as 
ot  ■*  00 ,  a  ramp  distribution  (beta  distribution  with  shaping  parameters  a  =  2, 
b  -  1  —  see  reference  1)  becomes  an  excellent  approximation  of  the  truncated 
Rayleigh  function. 

The  evasion  region  density  function  ultimately  must  be  combined  with  the 
TMA  uncertainty  region,  which  is  Gaussian  distributed,  to  obtain  online 
measure-of-effectiveness  (MOE)  algorithms.  Recent  investigations  have  shown 
that  this  is  very  difficult  to  achieve  using  the  beta  function  evasion 
regions.  However,  the  evasion  regions  modeled  here  are  Gaussian  distributed, 
and  it  is  well  known  that  a  Gaussian  evasion  region  combined  with  a  Gaussian 
uncertainty  region  yields  another  Gaussian  probability  region.  The  principal 
advantage  now  is  that  an  analytical  formulation  exists  to  perform  the  task  of 
developing  real-time  MOEs .  The  only  disadvantage  is  that  target  evasion 
course  is  limited  to  a  uniform  density  function,  and  target  evasion  speed  is 
limited  to  either  a  truncated  Rayleigh  or  a  ramp  density  function. 

Future  work  will  address  the  development  of  the  interactive  man-machine 
interface  necessary  to  combine  both  a  priori  knowledge  and  heuristic  aspects 
of  the  problem.  These  heuristic  aspects  (e.g.,  operator  selection  of  the 
evasion  course  and  speed  as  well  as  confidence  levels  for  these  selections) 
must  be  transformed  into  mathematical  models  that  develop  the  time-varying 
probability  density  functions. 
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